Cluster model of glass transition in simple liquids 
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On the basis of microscopic statistical mechanics of simple liquids the orientational interaction 
between clusters consisting of a particle and its nearest neighbors is estimated. It is shown that 
there are ranges of density and temperature where the interaction changes sign as a function of a 
radius of a cluster. The model of interacting cubic and icosahedral clusters is proposed and solved 
in mean-field replica symmetric approximation. It is shown that the glass order parameter grows 
smoothly upon cooling, the transition temperature being identified with the temperature of the 
replica symmetry breaking. It is shown that upon cooling a Lennard- Jones system, cubic clusters 
freeze first. The transition temperature for icosahedral clusters is about ten per cent lower. So the 
local structure of Lennard- Jones glass in the vicinity of glass transition should be most probably 
cubic. 
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Despite the growing interest to physical properties of 
liquids and glasses (see, for example, for the recent re- 
views lill^yj), nature of structural glass transition is still 
puzzling. While the experimental and phenomenological 
knowledge of non-ergodic amorphous phases has consid- 
erably improved in last time jlj , progress in first-principle 
statistical mechanical studies of structural glasses is 
much more slow. 

One of the promising ways to study the structural 
glasses is to explore the analogy between the phase tran- 
sitions in spin glasses (which are well understood now 
0,3) and the structural glass transition. In some extent 
this is motivated by the fact that the structure of the dy- 
namical equations for the appropriate correlation func- 
tions of super-cooled liquids and p-spin spin glass model 
(p > 3) are identical in the mean-field (MF) approxi- 
mation 0- ■ so this model may be used to describe, 
at least qualitatively, the properties of structural glasses. 
The main obstacle in this way is the absence in struc- 
tural glasses of quenched disorder (in contrast to the spin 
glasses). Spin glasses are microscopically quite different 
from liquids and thus seem not suitable for their descrip- 
tion. Furthermore, this approach does not include any 
information on the local structure of super-cooled liquids 
and glasses. In this paper we propose the scenario where 
the analog of quenched disorder in Lennard- Jones system 
appears in natural way allowing to apply the methods of 
spin glass theory to investigation of the structural glass 
transition. 

The structure of supercooled liquids and glasses is 
incompletely understood even for simplest systems. 
Frenkel Q has suggested a qualitative picture of local 
structure of a dense supercooled liquid: he supposed that 
in small volumes it has a crystal-like structure. On the 
other hand, as was argued by Frank an icosahedral 
cluster of 13 particles has a significantly lower energy 
than more obvious arrangements having the symmetry 



of FCC or HCP crystals. It was inferred from the com- 
puter simulations of a dense supercooled Lennard- Jones 
liquid 9] that the local symmetry of simple liquids is de- 
termined by defective or fragmented icosahedral building 
blocks that exhibit a five-fold symmetry. However, as 
was shown later, in a supercooled Lennard-Jones liquid, 
small crystalline FCC clusters nucleate following the tem- 
perature quench It seems that the local structure 
of a supercooled Lennard-Jones liquid is a result of com- 
petition between FCC and icosahedral local symmetries 

The concept of interacting clusters was used in phe- 
nomenological theories of bond orientational order 0, 
H^ l or "orientational melting" 0,^3- The microscopic 
approach to description of the bond-orientational order 
(or hexatic phase in two dimensions) was developed in 
[ill UtI Hif . In the present Letter we use this ap- 
proach to analyze the intercluster interaction and to in- 
troduce the model of glass transition in a Lennard-Jones 
liquid which elucidates the structural properties of the 
corresponding glass. 

Our starting point is the expression for the free en- 
ergy of the system as a functional of a pair distribution 
function <72(ri,ro) which has the form [ill llq: 

F/k B T = J drdr oP52 (r,r )ln [(A 3 P S 2 (r, r )) - l] - 

n n + l f 

~22 (n + l)\ J 5 »+i( r i"- r »+i)32(ri,r )--- 



(n+1)! 
• g 2 {r n +i, r )dri • ■ ■ dr n+1 dr 

J $(r - Y )pg 2 {r, r )drdr . 



(1) 



Here Sk+i{ri---rk+i) is the irreducible cluster sum of 
Mayer functions connecting (at least doubly) k + 1 parti- 
cles, p is the mean number density, <I>(r — ro) - interpar- 
ticle potential, A = /i/(27rmfc s T) 1 / 2 . 
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In an isotropic liquid pair distribution function 
£ 2 (r,r ) depends on |r-r | only: 5 2 (r,r ) =p(|r-r |), 
where g(r) is the radial distribution function. In the state 
with bond orientational order rotational symmetry of the 
pair distribution function is broken: 



52(r,r ) = g(\r - r |) + 8g(r,r ). 



(2) 



6g(r,ro) has the symmetry of the local environment of 
the particle at i"o and may be approximately written in 
the form [HO]!: 



^(r,r ) = /(n)<y(r,-|r-r |). 



(3) 



Here determines the direction of the vector r — ro , and 
r s is the size of the cluster, which is approximately equal 
to the diameter of the first coordination shell. Function 
/(fi) gives the probability of cluster orientation and may 
be expanded in a series in spherical harmonics: 



oo / 



f(r,Q)=J2 E fim(r)Y lm (Q). 



(4) 



1=0 m=-l 



In expansion Q only the terms corresponding to the lo- 
cal cluster symmetry should be retained. Coefficients fi m 
are the order parameters for the phase transition to the 
anisotropic phase. 

To estimate the intcrclustcr interaction, let us expand 
the free energy up to the second order in 8g(r,ro). 
Omitting the O-independent terms, one has: 

AF/k b T = -\j r(ri,r o y 2 ) x 

x 5g(Ti,r Q )5g(Ti,To)dT\dT2, (5) 

where [HGjl: 

r(r 1 ,r?,r 2 ) = ^-^-— f 5 fc+1 (n...r fc+1 ) x 
x f(l r 3-r?|)---3(|i-fe+i -r?|)dr 3 ---dr fc+ i = 

= E^Ti r '^ r ') E iWfii)i&(n 3 ). (6) 



l=-r 



The angles fii and f^2 determine the directions of the 
vectors r = ri — ro and r' = r 2 — ro. 

After substituting © and © in © we obtain: 

1 00 A 

AF(r s )/k B T = - 2 E ^7TT r '( r - r ^ X 

i 

^ m (ni)^(n 2 )/(ni)/(0 2 )rfdOid0 2 = 

!=— JTi 

1 oo Z 

= -oE J '(^) E i/^i 2 - ( ? ) 



Here J,(r s ) = ^ r %{r s ,r s ). 

Function AF(r s ) may be interpreted as the mean- field 
orientational interaction energy of the system of clusters 
having the size r s . To get the full energy of the system 
one should integrate Q over the probability of finding 
the cluster with the size r s which is given by the function 
r 2 g(r) in the vicinit y of the first maximum. 

As was shown in there is a simple approximation 
for r(ri,r5,r 2 ) which gives rather good results for the 
Lennard- Jones potential $z,j(r): 

r(n,r?,r 2 ) = p(exp(-$ LJ (|n - v 2 \)/k B T) - 1). (8) 

Using Eq. (JSJ one can obtain the estimation for Ji(r s ) 
as a function of r s . Fig^ represents J/(r s ) for I = 4 and 
6 along with r 2 g(r) in the vicinity of the first peak. It is 
seen that Ji(r s ) changes sign. From Eq. (JJJ) and Fig. Q] 
one can conclude that the frustration appears as a result 
of variations in the sizes of clusters due to local density 
fluctuations. Such kind of behaviour leads to frustration 
which is analogous to that in spin glasses. 
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FIG. 1: Ji(r s ) for / = 4 and 6 along with r^g(r a ) in the 
vicinity of the first peak as functions of r s for dimensionless 
temperature T = 1.0 and density p = 1.2. 

To study the transition in the system of interacting 
clusters we introduce the simple lattice Hamiltonian:: 



1 oo I 

H = -o E E4 E u^muun, 



(9) 



<i^j> l=0 m=— ! 



l=—m 



The functions Ui m {VLi) the lattice harmonics for the space 
groups corresponding to the cluster symmetry. Taking 
into account that < Ui m (Cli) >= fi m , one can see that 
in the MF approximation the energy calculated from Eq. 

coincides with the intercluster energy Q under ap- 
propriate choice of J' . 
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To simplify the problem we neglect in Hamiltonian © 
all the terms except for the ones corresponding to the 
unit representation of the space group. Furthermore, we 
consider only the cases I = 4 and I = 6 which repre- 
sent the cases of cubic and icosahedral symmetries. This 
Ising-like model may be called a "minimal" model. The 
Hamiltonian of the minimal model has the form: 



h = -learnt 



(10) 



Functions U = U(<p, 9) are the combinations of the spher- 
ical harmonics. We will consider separately symmetries 
of "simple" cube (I — 4, m = 0, ±4), cube (I = 6, m = 
0, ±4) and icosahedron (I = 6, m = 0, ±5) correspond- 
ingly 0,0]. The trace in this case is defined as follows: 
Tr(. . .) = dip J™ dcos(9)(. . .). For example, for / = 4 
one has: 



U=U(cp,9) 



(Y ii (<p,d) + Y 4i (-< P) e))} (ii) 




T/J 



The interactions J.^ are such that the MF approxima- 
tion gives exact solution (an infinite-range interaction). 
It is easily seen that in the minimal model l|10fl with- 
out disorder in the framework of the MF approximation 
there is a first order phase transition to the state with 
bond orientational order. From Fig. Q]it is clear that, as 
the first approximation, J l i} may be random interaction to Z = 6, m = 0, ±4 and icosahedron to I = 6, m = 0, ±5. 
with Gaussian probability distribution. 



FIG. 2: Order parameters m and g 1//2 as functions of 
T/J; Jo = 0. Dotted vertical curve shows the tempera- 
ture T A _ T which separates stable and unstable regions of 
the replica symmetric theory; the glass transition occurs at 
T A _ T . "Simple" cube corresponds to / = 4, m = 0, ±4, cube 



P(Jij) = 



exp 



(Jij - J ) 



- 1 



2J 2 



(12) 



Ijl3p give the following equations for these order param- 
eters: 



with J = J/VN , J = J /N. 

Then the free energy in the replica-symmetric approx- 
imation is equal to [lSI . |20| : 



(CO, P={U 2 ), q={Uf 



(14) 



F 




J0_\ ^ +t 2<f_ 

kT 2 4 



exp 



InTr 



4 

exp 



(') 



(13) 



Here t = J/ksT and 



ztyjq- 



Jo_ 

kT 



where (...) = Tr(. . . e e )/Tre e and (...) = 
j^ 3 ^ -y=e~ z / 2 [. . .]. The temperature dependence 
of the order parameters is represented in Fig. [2] One can 
see that the glass phase grows smoothly upon cooling. 
The similar behavior takes place in a quadrupolar glass 
[3113,111. It is not clear now whether the appearance 
of the nonzero bond-orientational order is an artifact of 
the model or an inherent property of the glass transition 
in the Lennard- Jones system. 

We define the transition temperature as the replica 
symmetry breaking temperature. The replica symmetric 
solution is stable unless the replicon mode energy A rep i is 
nonzero |2Ct |22j . For our model we have 



The order parameters are: m is the regular order pa- 
rameter (an analog of magnetic moment in spin glasses), 
q is the glass order parameter and p is an auxiliary order 
parameter. The extremum conditions for the free energy 



A rep l = 1-* 2 ((C/ 2 )) 5 



(15) 



where ((■■■)) denotes the irreducible correlator. The tem- 
perature T A _ T corresponding to A rcp i = defines the 
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glass transition. In the Fig. [3] we plot the dependence 
T A _ T on J for I = 6. From Fig. |3] one can conclude 
that in the vicinity of the glass transition the cubic lo- 
cal symmetry appears first and, taking into account the 
contribution of the terms with I = 4, is preferable for 
the Lennard- Jones glass. This conclusion is in agreement 
with the result in Ref. 0] > where it was shown that the 
correlation length of the cubic bond orientational order 
exceeds the icosahedral one. 

From the condition of marginal stability 0, E| we also 
evaluated numerically the so-called dynamical transition 
temperature To at Jo = and found that within the ac- 
curacy of calculations Td and T A _ T coincide. One can ex- 
pect that in the minimal model (jlOJl full replica symmetry 
breaking should occur (as in the Sherrington-Kirkpatrick 
model jj|) but further investigation is necessary to eluci- 
date this question. 
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FIG. 3: Evolution of the temperature T A _ T with Jo /J 
is shown in the figure for the symmetries of cube and 
icosahedron for I — 6; below T A T the replica sym- 
metric MF solution becomes unstable. The free energy 
F(T A _ T (Jo/ J))/NkBT A _ T is shown in the inset. 

In summary, using the approach developed previously 
for study of a bond-orientational order in simple liquids 
we estimated the orientational interaction between 
clusters consisting of a particle and its nearest neighbors. 
It is shown that there are ranges of density and temper- 
ature where the interaction changes sign as a function of 
a radius of a cluster, the probability of a radius being 
determined by the radial distribution function of a liquid 
in the vicinity of the first peak. The model of interacting 
cubic and icosahedral clusters was solved in mean-field 
replica symmetric approximation. Due to absence of re- 
flection symmetry 20] the glass order parameter grows 
smoothly upon cooling. The transition temperature is 
identified with the temperature replica symmetry break- 



ing. It was shown that upon cooling a Lennard- Jones 
system cubic clusters freeze first. The transition tem- 
perature for icosahedral clusters is about ten per cent 
lower. So the local structure of Lennard- Jones glass in 
the vicinity of glass transition should be most probably 
cubic. 
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